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Abstract. I present a full leading-order calculation of F2{x, and Fl{x, Q^), 
including contributions not only from leading order in as , but also from the lead¬ 
ing power of Os for each order in ln(l/x). The calculation is ordered according 
to the inputs and evolution of the structure functions, and the perturbative form 
of the inputs is determined. I compare the results of fits to data to those using 
conventional LO and NLO order calculations, and the correct inclusion of leading 
ln(l/a;) terms is clearly preferred. A prediction for Fl{x, is produced which 
is smaller at small x than that obtained from the conventional approach. 

INTRODUCTION 

There has recently been a great deal of new data obtained at HERA for the 
structure function F 2 {x,Q‘^) at small x [1,2], and consequently a great deal 
of theoretical activity. The main theoretical question is whether one should 
include leading ln(l/a;) terms which cause small x enhancement of terms at 
high orders in ag, and if so, then in precisely what manner. Common wisdom 
was that such terms should not be included because standard ways of doing 
so showed that the terms were indeed important for the ranges of x and 
probed at HERA, but that they worsened global fits to data obtained from 
conventional LO and NLO approaches, rather than improved them. Indeed, 
the best global fits to data seemed to come from NLO evolution starting at 
2GeV^ where the input for the singlet quark distribution (weighted by 
x), and hence the structure function, behaved like at small x [3]. Thus, 

despite the good quality of the global fits, the situation was unsatisfactory 
since terms were ignored because they were inconvenient, rather than for any 
sound theoretical reason, and the required input for the quark distribution 
was of a rather steep, unjustified form. 
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Recently I proposed that a correct theoretical treatment of the calculation 
of structure functions demands that the leading powers of as for given powers 
of ln(l/a;) must be included, but in expressions for physical quantities [4], This 
latter point was inspired by, and is similar to the idea of using physical anoma¬ 
lous dimensions proposed by Catani [5]. However, the structure imposed by 
an ordered expansion within a given renormalization scheme has a number 
of other consequences for the form of the expressions for structure functions. 
I will discuss the leading-order-renormalization-scheme-consistent (LORSC) 
calculation of structure functions briefly before looking at comparison with 
experiment. 


THE LORSC CALCULATION. 

The principles underlying the calculation and the consequences are: 

1. The quantities one calculates to a given order are directly observable. 
Hence rather than coefficient functions and parton distributions the structure 
function is factorized into inputs at some scale Qj and evolution away from this 
scale. One obvious consequence of this is factorization scheme independence. 

2. For a physical quantity the leading-order expression for each independent 
component begins at its lowest power of a*, whatever power that happens 
to be, e.g. if a term of ln™'(l/a;) appears for the first time at order a” this 
is the leading-order term for x-dependence of this type. In practice this is 
implemented in moment space where the leading power in for each inverse 
power of the moment variable N is part of the full leading-order expression. 
As it must, this results in leading-order expressions which are renormalization 
scheme independent and which are therefore compatible with the use of the 
one-loop running coupling. 

3. The inputs for the structure functions are two fundamental nonperturbative 

functions, one each for F 2 and Fl, convoluted with perturbative contributions. 
The nonperturbative functions are taken to be flat at small x. The perturba¬ 
tive parts are determined so that the total expression for the structure function 
is independent of the choice of starting scale, i.e. the LORSC expression is 
invariant under changes of Qj up to changes beyond leading order. This deter¬ 
mines the inputs uniquely up to a scale, Q%p, which should be typical of the 
onset of nonperturbative physics. For both F 2 and Fl the inputs are required 
to behave roughly like for Q/ ~ 20 — lOOGeV^ for 10“^ < a; < 10“®. 

Moreover, the forms of Fl{Q‘])-, F 2 {Q‘]) and {dF 2 {Q^)/d \nQ‘^)Q 2 are all re¬ 
lated at small x, and once one of them is set (in practice F 2 {Qj)) by htting to 
data there is very little freedom in the others, a constraint largely absent in 
conventional approaches. 

These points are discussed in far greater detail in [4], and the way in which 
they are put into practice is presented in the latter of these papers. 
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FITS TO DATA AND PREDICTIONS 

The fit using the LORSC expressions is performed with the parameters 
specifying the form of the nonperturbative inputs and the nonperturbative 
scale left free, and is repeated for a wide variety of input scales Q/- The best 
£t comes from = 40GeV^, but is very insensitive to Q] for values from 
20GeV^ ^ 80GeV^. We find that = O.SSGeV^, precisely the sort of 
valne expected. The charm and bottom quarks are dealt with rather naively, 
treated as massless above threshold (4GeV^ and 20GeV^ respectively) and 
playing no role below this. The charm threshold is chosen to give a reasonable 
description of data on the charm structure function. is chosen as the 
(squared) renormalization scale, and the value of Aqcd is held fixed at lOOMeV 
(giving as{Mz) = 0.115 at one loop), though this is certainly close to the best 
fit value. The results of the fit are shown in table 1. 

The conventional NLO fit is performed in the usual manner and the input 
scale for the parton distributions Qq is chosen to be the same as that for the 
charm threshold (2.75GeV^). For the best fit we find that A|^ = SOOMeV, i.e. 
as{Mz) = 0.118. The input quark distribution ~ The conventional 

LO fit is performed in the same way with charm threshold of 3GeV^, Aqcd = 
lOOMeV, (one-loop as{M‘^) = 0.124), and small x qnark ~ The 

results of these more conventional fits are also shown in table 1.^ 

In each case the values of F 2 {x, Q"^) used for the HERA data are not simply 
those published but are corrected for the values of Fp(a;, Q"^) predicted in each 
of the fits. For the LORSG fit this causes a slight lowering of F 2 {x,Q'^) for 
some points due to smaller predicted Fl{x, Q^) which in fact improves the fit 
to the small x data by a of abont 5. For the NLO fit there is essentially no 
correction, but for the LO fit the predicted Fl{x, Q^) is extremely large since 
the fit requires both a large coupling and large gluon at small x. Hence, the 
values of F 2 {x,Q'^) are corrected upwards, sometimes by 4%, leading to the 
correct LO fit being ~ 25 worse than a fit to nncorrected data.^ 

As one can see the LORSG fit is clearly superior: although the difference 


All fits were constrained not only by consistency with charm structure function data, 
but also by consistency with high-x prompt photon data. 

For the LO and NLO fits the normalization of the HI data is as low as I allow, i.e. 
98.5%, while for the LORSC fit it is 100%. In each case the normalization of the ZEUS 
data is 1% higher. 


TABLE 1. Quality of fits using full leading-order (LORSC) 


expressions and conventional LO and NLO expressions. 


x-range 

data points 

LORSC 

NLO x^ 

LO x" 

X > 0.1 

551 

622 

615 

605 

X < 0.1 

548 

483 

554 

598 

total 

1099 

1105 

1169 

1203 
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in quality is ~ 0.06 per point this is significant when considering 1099 data 
points. Indeed, the difference between the quality of the LORSC £t and the 
NLO fit is larger than that between the NLO and LO fits. The quality in 
terms of different experiments may be found in the second of [4] (as may the 
references for the experiments), but it is illustrative to present it in terms of 
two bins in x. There is little difference in the quality of the high-x fits, as we 
would expect since at high x the missing ln(l/a;) terms in the LO and NLO 
fits cause no problems (and an appropriate choice of coupling makes the LO 
fit as good as the NLO fit). In fact, the LORSC calculation is missing NLO 
as terms, and should perhaps be a little worse than the NLO ht.^ However, 
at small x there is a very clear deterioration in the quality of fit going from 
LORSC to NLO to LO (which would be even greater if values of A consistent 
with the best high-x hts were used). This is due to the loss of vital ln(l/a;) 
terms as ones goes from LORSC to NLO, and also, to a lesser extent, when 
going from NLO to LO. Hence, the results of the fits are precisely as we would 
expect from theoretical arguments. They clearly imply that in determining 
values of q;s(M|) and the gluon distribution in particular schemes the NLO 
fit should only be trusted at relatively high x. 

Shown in fig. 1 is the comparison of the predictions for Fl{x, Q^) made us¬ 
ing the best fit from the LORSC approach and the NLO approach. As one can 
see there is significant difference between them, and any relatively accurate 
measurement of Fl{x, should see some sign of which approach is preferred. 
The “determination” of Fl{x,Q‘^) by HI is no use in differentiating between 
the two since it is really just a consistency check. If one assumes that the HI 
NLO fit is correct and extrapolates from low y = jxs, where the determi¬ 
nation of F 2 {x, Q'^) from the cross-section is insensitive to Fl{x, Q^) then the 
difference between the measured cross-section and the extrapolation provides 
a value of Fl{x, Q"^) [6]. This must be consistent with the NLO prediction for 
Fl{x, Q‘^), which it is. However, assuming that the LORSC calculation is cor¬ 
rect the extrapolation is different, and hence so is the determined Fl{x,Q‘^). 
Consistency between the determined Fl{x, and the prediction is also good 
in this case. ^ Hence, a true measurement of Fl{x, Q"^) would be an invaluable 
aid to the determination of the real physics describing hadron interactions at 
small X. An important role may also be played by improved data on the charm 
structure function, and incorporating heavy quarks correctly into the LORSC 
approach is a project well underway. Less inclusive quantities may well also 
play an important part, but calculations remain to be done. 

This is clear if we perform fits to high-x data only: the LORSC fit only improves by a 
few points, but the NLO fit improves to 578, with = 215MeV. Such a low value of 
is not allowed in the full NLO fit since it leads to a very poor ht at small x. 

The HI fit results in a gluon which is inconsistent with both high-a; prompt photon data 
and moderate-a: charm data, and it would be illustrative to repeat the process with a more 
constrained NLO ht. It would also be interesting to use all available high y data, i.e. also 
include ZEUS data. Such a study is being performed. 
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FIGURE 1. Comparison of predictions of Fl{x,Q^) using the LORSC fit and the NLO 
fit. In both cases Fl{x, Q^) increases with at small x. 


CONCLUSION 

A comparison to all significant data gives clear evidence that the LORSC 
calculation of structure functions is preferred to the usual NLO-in-a^ ap¬ 
proach. Thus, although the usual expansion technique seems to be acceptable 
at small x at present, and will probably be used for most QCD calculations 
in the near future, when attempting to describe any data at small x it should 
be borne in mind that this approach may well be untrustworthy: anomalies 
may occur, and more careful calculations may be needed. The situation will 
be clarihed by a variety of different small x measurements. 
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